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Anomalous Field- Angle Dependence of the Specific Heat of Heavy-Fermion Superconductor UPts 
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We have investigated the field-angle variation of the specific heat C(H, (f>, 6) of the heavy-fermion su- 
perconductor UPt3 at low temperatures T down to 50 mK, where ip and 9 denote the azimuthal and polar 
angles of the magnetic field H, respectively. For T = 88 mK, C(H,6 = 90°) increases proportionally to 
yfn up to nearly the upper critical field /fc2, indicating the presence of line nodes. By contrast, C{H, 9 = 0°) 
deviates upward from the ^fH dependence for yftTjll^ > 0.5. This behavior can be related to the suppres- 
sion of Hc2 along the c direction, whose origin has not been resolved yet. Our data show that the unusual 
H^2 limit becomes marked only when 9 is smaller than 30°. In order to explore the possible vertical line 
nodes in the gap structure, we measured the (f> dependence of C in wide T and H ranges. However, we did 
not observe any in-plane angular oscillation of C within the accuracy of AC/C ~ 0.5%. This result implies 
that field-induced excitations of the heavy quasiparticles occur isotropically with respect to 0, which is 
apparently contrary to the recent finding of a twofold thermal-conductivity k{^) oscillation. 
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1. Introduction 

UPt3 is an unconventional superconductor with multiple 
superconducting (SC) phases [1) (A, B, and C phases, see 
Fig. 1). In particular, strong evidence for spin-triplet super- 
conductivity has been provided by NMR Knight shift mea- 
surements, [2, 3) and the presence of line nodes in the SC 
gap has been clearly indicated by the power-law behavior in 
thermodynamic and transport quantities. [4-8) To explain the 
multiple SC phases, the presence of a two-component order 
parameter coupled to a symmetry-breaking field (SBF) has 
been proposed; whereas two order parameters are degener- 
ate in the B phase, one of them appears in the A phase and 
the other appears in the C phase. Although the SBF, which 
lifts the degeneracy of the order parameters at zero field, is 
still controversial, it could possibly be related to short-range 
antiferromagnetic (AFM) fluctuations in the basal plane de- 
veloping below 5 K [9-11) and/or distortions in the crystal 
structure [12-14). 

While UPtj has been well established to have nodes 
in the SC gap, the precise gap structure of each SC phase 
has been debated actively. On theoretical grounds, various 
multi-component order parameters have been proposed so 
far. [1, 15, 16) One of the plausible scenarios is the E2u model 
with spin-triplet pairing involving the SC order parameters 
- A:^) and Ik^kyk^, whose gap functions have fourfold 
symmetry with vertical line nodes and are degenerate in the 
B phase as \k;\{k\ + A;^) with axial symmetry. Experimen- 
tal evidence supporting this scenario has been provided by 
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small-angle neutron scattering [17, 18) and junction experi- 
ments [19, 20) performed in the A and B phases. However, 
a conflicting scenario, the spin-triplet ^i,, model with the or- 
der parameters k^{5k^ - 1) and k^iSkl - 1) has recently been 
proposed on the basis of the results of field-angle-resolved 
thermal-conductivity K{(p) measurements performed in the B 
and C phases; [2 1 , 22) under a magnetic field H rotated around 
the c axis, k{^) exhibits a twofold oscillation at 50 mK only in 
the C phase. 

In addition to the SC gap structure, the origin of the 
strong limit of the upper critical field //c2 for H II c at low 
T has not been resolved. Although it resembles the //c2 limit 
due to the Pauli paramagnetic effect, [23, 24) Hd cannot be 
affected by this effect in UPt3 because the results of NMR 
Knight shift experiments suggested that the d vector of the 
triplet Cooper pair is perpendicular to H in the C phase; the 
Cooper pairs can be polarized along the field direction. [2) 
This is reminiscent of a similar, unusual H^^o limit observed 
in another plausible spin-triplet superconductor Sr2Ru04 for 
H\\ab. [25,26) 

To obtain hints for clarifying the precise gap structure and 
the origin of the H^^i limit in UPt3, we have performed field- 
angle-resolved specific-heat C{H, (p, ff) measurements down 
to 50 mK, where (f> is the azimuthal field angle within the ab 
plane measured from the crystalline [1210] axis and 9 denotes 
the polar angle between H and the c axis. Because C/T is pro- 
portional to the quasiparticle (QP) density of states (DOS) at 
low T, this technique enables us to get an insight into nodal 
directions of bulk superconductors from the field-orientation 
dependence of QP excitations. [27) At 88 mK, an unusual 
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Fig. 1. (Color Online) Schematic field-temperature phase diagram of UPtj 
for H II ab. Squares (circles) represent the field and temperature at which 
C(0) measurements are performed with 6 = 90° (60°). The inset represents 
the 9 dependence of H^i determined from the C(H) data obtained at 88 and 
150 mK. Dashed lines are guides to the eye. 



increase in C{H) for -\/hJH^ > 0.5 was found only when 
\6\ < 30°. Correspondingly, the dependence of C at 88 mK 
exhibits a peak at = 0° with a width of 16*1 < 30° in the 
high-// regime, whereas C{0) at low H shows a twofold oscil- 
lation with a deep minimum at = 0°. This feature indicates 
the occurrence of unusual QP excitations for l^l < 30° in the 
high-// regime, which might be related to the //c2 limit. Un- 
der H rotated around the c axis, we observed no variation in 
C{(p) down to 50 mK (~ OAT^), implying a small change in 
QP DOS as a function of in-plane field angle. 

2. Experimental Details 

A single crystal of UPt3 used in the present study (Tc = 
0.52 K, 54.2 mg) was grown by the Czochralski pulling 
method. The sample was shaped into a disk with flat surfaces 
parallel to the ab plane. The specific heat C was measured by 
relaxation and standard quasi-adiabatic heat-pulse methods in 
a dilution refrigerator (Oxford Kelvinox AST Minisorb). In 
all the data presented, a nuclear Schottky contribution (Cn) 
by the Zeeman spUtting of '^^Pt nuclei (/ = 1/2) and the ad- 
denda contribution were subtracted. It was confirmed that the 
addenda specific heat is always less than 4% of the sample 
specific heat and exhibits a negligible field-angle dependence. 
Although a rapid increase was observed in C/T at zero field 
below 0.1 K, [28,29) it is probably related to an AFM long- 
range ordering at 20 mK [30,31). Magnetic fields were gener- 
ated in the xz plane using a vector magnet composed of hor- 
izontal split-pair (5 T) and vertical solenoidal (3 T) coils. By 
rotating the refrigerator around the z axis using a stepper mo- 
tor mounted at the top of a magnet Dewar, we controlled the 
field direction three-dimensionally. Thus, we performed ex- 
periments in an accurate (< 0.1°) and precise (< 0.01°) field 
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Fig. 2. (Color OnHne) (a) Magnetic-field dependence of (C - Cn)/r at 
several polar angle S values at 88 mK. (b) (C - Cn)/r as a function 
of ^HIH^2(8). Dashed lines ai'e guides to the eye. Solid lines represent 
the fits to the data using the function a Vff + h in the low-field regime 
iiioH < 0.5 T). 



alignment with respect to the ab plane, which was adjusted 
by using the anisotropy of the specific heat of UPt3 in the ac 
plane. For each point, the specific-heat value was determined 
by the average of about ten successive measurements. 

3. Results and Discussion 

3. 1 Field dependence of the specific heat 

Figure 2(a) shows the H dependence of the nuclear- 
subtracted C divided by T at various 6 values. These data were 
obtained at 88 mK and = 0°. At low H below 0.5 T, C{H) 
increases rapidly with H at any 9. This rapid increase is fit- 
ted well using the function a y/H + b, as represented by solid 
lines in Fig. 2(b), where C/T is plotted as a function of the 
square root of H/Hc2(6)- The occurrence of the V// behavior 
at any indicates that Doppler-shifted QP excitations around 
line nodes are predominant in determining C(H) and excludes 
the possibility of forming a SC gap structure with only point 
nodes at the poles. [32) 

If //c2 is dominated by the orbital pair-breaking effect, the 
C{H) of anisotropic superconductors with line nodes is ex- 
pected to exhibit the V^-linear dependence up to nearly Hc2- 
[33) Such behavior is indeed observed in C{H) at 6 - 90°. 
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By contrast, at 6 - 0°, C as a function of exhibits a 
gradual upturn in the high-// regime for ^ 0.5, as 

shown in Fig. 2(b). This upturn is clearly seen for \6\ < 30° 
and enhances with decreasing 6 to zero. Although the y/H- 
linear dependence at low H is slightly smeared at any 9 by 
the effect of thermal QP excitations, a similar feature was ob- 
served at 150 mK as well (not shown). The characteristic field 
of 0.25//c2, above which the upturn becomes marked, is much 
lower than the transition field from the B phase to the C phase 
of about 0.7//c2 for H || c [34). 

The field-angle 9 dependence of Hc2 determined from the 
C{H) data obtained at 88 and 150 mK is shown in the inset 
of Fig. 1. Also, a relatively strong Hc2 limit occurs for \9\ < 
30°. Thus, the anomalous H dependence of C{H) observed 
for \9\ < 30° seems to be related to the unresolved mechanism 
of the //c2 limit. 

A possible origin of the H^i limit is the Pauli paramag- 
netic eff'ect; the H dependence of C at = 0° is qualitatively 
similar to the calculated result based on a quasiclassical the- 
ory with the paramagnetic parameter ju - 0.86 for nodal d- 
wave superconductors. [33) In addition, the multi-band effect 
cannot be ruled out because UPt3 has five bands in the Fermi 
surface [35); a possible change of the relevant band that deter- 
mines //c2 when T, H, and field orientation are varied could 
lead to the apparent //c2 limit. Nevertheless, no strong evi- 
dence has yet been established for both possibilities. In partic- 
ular, the Pauli effect is incompatible with the results of NMR 
Knight-shift experiments. [2, 36) Thus, the origin of the Hd 
limit remains unclear. 

3.2 Field-angle 9 dependence of the specific heat 

The 9 dependence of C measured under H rotated within 
the ac plane {(p - 0°) is shown in Fig. 3, where the C{9) data 
obtained at the interval 0° < 6* < 90° are plotted repeatedly 
in accordance with the crystal symmetry. Let us first focus 
on C{9) in the normal state. At 600 mK, which is sufficiently 
above T^, C is invariant under the 9 rotation at any H below 
2.8 T. By contrast, at low T below 300 mK, C{9) at 2.8 T 
that is higher than Hd exhibits a small twofold oscillation 
with a maximum at 9 - 90° (H \\ a). This anisotropy might 
originate from the anisotropy of the DOS enhancement by a 
metamagnetic transition occurring at 20 T for H \\ ab. [37) 

In the SC state, C varies markedly with changing 9. 
At low T below 150 mK and low H below 1 T, C(9) ex- 
hibits a large twofold oscillation with a deep minimum at 
9 ^ 0° (H II c). The relation C(9 = 0°) < C(90°) means 
smaller QP excitations at low H in H \\ c rather than in 
H II a. This anisotropy is opposite to the expectation from 
the Hc2 anisotropy {Hd\\a > Hc2\\c), which naturally leads to 
C(0°) > C(90°). This result implies that, if the limit is 
absent, the intrinsic Hd(T = K) is larger for H \\ ab than 
for H \\ c, the same anisotropy of observed in the vicinity 
ofTe. [38). 

By increasing H at 88 mK, a small peak whose width is 
|6i| < 30° develops in C{9) at 6» = 0°. At 1.8 T, this peak makes 
the anisotropy in C{9) reversed [C(0°) > C(90°)], consistent 
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Fig. 3. (Color Online) Polar-field-angle dependence of (C-Ca)/T at sev- 
eral temperatures and fields. Dashed lines are guides to the eye. For easy 
viewing, the data measured at the interval 0° < d < 90° are plotted repeat- 
edly (6 is converted to -6 and 180° ± 6). 



with the expectation from the observed anisotropy. A sim- 
ilar peak can also be seen at 150 and 300 mK. In particular, at 
300 mK, the peak is visible even at low H and with a larger 
width (|6I| < 45°). The observed peak in C(9) at 88 mK corre- 
sponds to the anomalous increase in C{H) for \9\ < 30° shown 
in Fig. 2(b). Therefore, this peak might define the conditions 
of the unresolved mechanism inducing the unusual Hc2 limit 

of UPt3. 

3.3 Field-angle (p dependence of the specific heat 

We next focus on the variation in C under H rotated 
around the c axis, C{<p). This C{4>) measurement is a powerful 
tool for identifying the SC gap structure. According to the re- 
sults of the Doppler-shift analysis assuming an anisotropic SC 
gap on a quasi-two-dimensional Fermi surface, zero-energy 
DOS, which corresponds to C/T{cf>) at low T, is expected to 
be enhanced (suppressed) when H is applied along the antin- 
odal (nodal) direction. [39^1) Even at a moderate temper- 
ature, a clear C(0) oscillation reflecting the gap structure is 
expected to be observed although its sign changes with the 
vortex-scattering effect. [42,43) Indeed, a fourfold oscillation 
reflecting the iii2_,,2-wave gap has been observed in the C(<^) 
of CeCoIns [44) and Celrlns [45) in the wide T and H regions. 

Contrary to the expectation, however, no angular varia- 
tion was detected within an experimental error in the C(0) 
of UPt3 under an in-plane rotation of H (9 - 90°) at tem- 
peratures 50 < r < 200 mK. The T-H ranges in which we 
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Fig. 4. (Color Online) Azimuthal-field-angle <f> dependence of (C - C„)/T. 
Bars represent 1% of Ch/T = [C(H) - C„(H) - C(0)]/r. 



searched for the C{(p) oscillation at = 90° are summarized 
in Fig. 1 with squares. For instance, the C(0) data measured 
in the normal state at 88 mK and 2.8 T, and that measured in 
the SC state at 50 mK and 1.5 T are shown in Figs. 4(a) and 
4(b), respectively. The C(4>) data measured under other T-H 
conditions are also presented elsewhere. [29) 

On the basis of our fits to the data using the function 
Co(r) + Ch{T,H){1 - A„cos«0) (« = 2 or 4), where Co 
and Ch are the zero-field and field-dependent components of 
the electronic specific heat, respectively, and A4 (A2) is the 
amplitude of the fourfold (twofold) oscillation normalized by 
Ch, A„ never becomes larger than 0.3% in the present T and 
H regimes. This is in sharp contrast to the observation of a 
twofold oscillation in the /c(0) of UPt3 at 50 mK, whose am- 
plitude is about 3% with respect to the normal-state value 
K„ [21) as well as a fourfold oscillation in C((p) with A4 larger 
than 3% at 0.4rc in the c/j2_,.2-wave superconductor Celrlns 
(re = 0.4 K). [45) 

If the SC gap has a horizontal line node at A:- = (e.g., 
the E2u model), the anisotropy of C(0) could be suppressed 
significantly at - 90° because H is always in the nodal di- 
rection for - 90° [see Fig. 4(e) of ref. 45]. To examine this 
possibility, we measured C{(p) under a conically rotating H ti- 
tled by 60° from the c axis (6 = 60°). The T and H conditions 
under which we investigated C(0) at 6* = 60° are denoted in 
Fig. 1 by circles. Nevertheless, as presented in Fig. 4(c), no 
angular variation was detected in C(0, 9 = 60°) at 88 mK in 
0.5 < fi()H < 2.2 T. These findings lead to a conclusion that 
QP excitations detected from the C{(p) measurements hardly 
depend on the field orientation around the c axis. We also ob- 
served that C{4>) measured after zero-field cooling and that 
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after field cooling in 4 T applied above 10 K along one of the 
three equivalent a axes provided the same results. 

For nodal superconductors, QP DOS is expected to 
be scaled using a function of a single parameter x = 
(T/Tc)/ -\/H/Hc2 [46-48). Theoretical calculations predicted 
the scaling relations of various physical quantities, such as 
C(H,T)/T^ oc f(x) and k(H,T)/T^ oc g(x) at low T, where 
f{x) and g{x) are constant for a large x and proportional to l/x 
for a small x. Indeed, C(H, 0) oc ^/7^/^^^ and C(0, T) oc T-, 
which are expected for nodal superconductors, can be given 
by this scaling relation. For UPt3, it was reported that k/T^-"^ 
in H II ab can be scaled by x at low H and low T. [49) We 
found that C{T, H) divided by r'^- can also be scaled by x in 
the intermediate-T regime. Figure 5 shows C/T^-^^ as a func- 
tion of l/x plotted using the C{T) and C{H) data at 6* = 90°. 
In 150 < r < 220 mK, C/r'^^ proportional to 1 /x when x 
is small and tends to saturate when x is large. These features 
support that the specific heat of UPts reflects QP DOS excited 
around the nodes. 

However, at lower T below 88 mK, C/T^'^^ does not fol- 
low the scaling relation owing to the unusual increase in C{T) 
at low T . Although one might consider that the variation in 
C by H, (p, and 9 is dominated by this unusual phenomenon, 
it should be noted that C{H) at 60 and 88 mK still shows the 
V^-linear dependence at low H (see Fig. 5), and C in the nor- 
mal state is known to be robust to the intensity of // up to 7 T 
down to 50 mK. [50,51) These findings imply that QP excita- 
tions still play an important role in determining the change in 
C{H, (f>, 9) at low T. 

In order to settle the issue of the SC gap structure of 
UPt3, the apparent discrepancy between the absence of any 
C(0) oscillation and the presence of a twofold oscillation in 
Kief)) should be resolved. A similar discrepancy has been re- 
ported on a filled skutterudite superconductor PrOs4Sbi2 as 
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well. [27, 52) The discrepancies might be explained by con- 
sidering the multi-band effect because C{H, (p, 9) mainly de- 
tects the contribution of a heavy-mass band whereas k{H, (p, 9) 
may not. Indeed, UPt3 has four heavy bands and one relatively 
light band. [35) In this context, the results of both C{(p) and 
K{(f)) experiments become consistent with each other if verti- 
cal line nodes exist only in the light band. Such assumptions, 
however, seem to be very speculative at present because of 
the lack of experimental and theoretical justifications. Further 
efibrts are needed to solve this issue. 

4. Summary 

We have measured the specific heat of UPt3 at low tem- 
peratures down to 50 mK in various field orientations. Un- 
der a magnetic field rotated in the ac plane, the specific heat 
exhibits a large twofold oscillation at low temperatures and 
low fields with a deep minimum at 9 - 0° (H \\ c), possibly 
reflecting the anisotropy of the intrinsic orbital limit of Hd- 
With increasing field or temperature, a striking peak develops 
at 61 = 0° with a width of \9\ < 30°. Correspondingly, at low 
T, an abnormal increase was observed in C{H) for l^l < 30°. 
Although its origin remains unclear, it might give the con- 
ditions in which the unusual H^o limit becomes prominent 
in UPta. Under a magnetic field rotated around the c axis in 
50 < r < 200 mK, the specific heat hardly depends on the 
field orientation. Further experimental and theoretical inves- 
tigations are essential to settle the SC gap structure of UPt3 
as well as the origin of the anomalous phenomena found for 
\9\ < 30°. 
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